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Abstract 

We present an updated analysis of the rare decay K L — > -n°ji + n~ within the Standard 
Model. In particular, we present the first complete calculation of the two-photon CP- 
conserving amplitude within Chiral Perturbation Theory, at the lowest non-trivial order. 
Our results confirm previous findings that the CP-conserving contribution to the decay 
rate cannot be neglected. By means of an explicit two-loop calculation, we show that 
this contribution can be estimated with sufficient accuracy compared to the CP-violating 
terms. We predict B(K L — > 7r /! + /! - ) SM = (1.5 ± 0.3) x 10~ n , with approximately equal 
contributions from the CP-conserving component, the indirect-CP-violating term, and 
the interesting direct-CP-violating amplitude. The error of this prediction is mainly of 
parametric nature and could be substantially reduced with better data on the K s — > 
7r°£ + £~ modes. The Standard Model predictions for various differential distributions and 
the sensitivity to possible new-physics effects are also briefly discussed. 



1 Introduction 

The rare decays Kl — > tt°£ + £~ are among the most interesting channels for precision 
studies of CP violation and flavour mixing in the AS = 1 sector p. The decay am- 
plitudes of these processes have three main ingredients: i. a clean direct-CP-violating 
component determined by short- distance dynamics; ii. an indirect-CP-violating (CPV) 
term due to K^-Kg mixing; in. a long-distance CP-conserving (CPC) component due to 
two-photon intermediate states. Although generated by very different dynamics, these 
three components are of comparable size within the Standard Model (SM). The precise 
knowledge about their magnitude (particularly of ii. and Hi.) has substantially improved 
in the recent past [2 . This improvement has been made possible by the observation of the 
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Ks — > vr°£ + £~ decays [3J Ej and also by precise experimental studies of the Kl — > 7r°77 
diphoton spectrum [SJ E] • 

The main difference of the — > 7c°fi + fi~ mode with respect to the electron channel, 
which has recently been re-analysed in Ref. |2j, is the two-photon CPC contribution. The 
latter corresponds to transitions into final states where the lepton pair has J CP = ++ 
or 2 ++ . The J CP = ++ case can safely be neglected in the electron mode, because of the 
strong helicity suppression pJIHJEl- On the contrary, the K L — > tt°(77) j=0 — ► n°(£ + £~) J=0 
amplitude generates the by far dominant CPC contribution in the muon mode [El EH- The 
main purpose of the present paper is a detailed analysis of this CPC amplitude, which is 
a fundamental ingredient to obtain a reliable SM prediction for the entire Kl — > 7r°// + /i~ 
decay width. 

Chiral Perturbation Theory (CHPT) provides a natural framework to analyse the CP- 
conserving Kl — > 7r°/i + /i~ amplitude. Here the first non-trivial contribution arises at 
0(p 4 ) from two-loop diagrams of the type Kl — > tt°(P + P~) — > 7r°(77) — > 7T°fi + fi~, where 
P = 7i, K. This structure is very similar to the two-loop Ks — * [i + amplitude analysed 
by Ecker and Pich jllj: there are no counterterms at 0(p A ) and the sum of all loop 
diagrams yields a finite unambiguous result. The Ks — > /i + 'pT amplitude of Ref. [TT] has 
indeed been used by Heiliger and Sehgal JUj to estimate the dominant pion-loop terms in 
Kl — > vr°/i + yU~. These authors extracted from the Ks — > ^ fi~ amplitude an appropriate 
electromagnetic form factor, describing the (7r + 7r~)j = o — > (77) — > transition, and 

used it in conjunction with a constant — > 37r vertex to estimate the Kl — > 7r° / u + /i~ 
amplitude. Although quite reasonable from a phenomenological point of view, this result 
is not fully justified within CHPT. 

By means of an explicit two-loop calculation, we show that the factorization of the 
electromagnetic form factor holds exactly at 0(p 4 ), independently of the structure of the 
K L — > 3n amplitude, and independently of the nature of the charged pseudoscalar meson 
inside the loop. We are then able to present the first complete 0(p 4 ) estimate of the CPC 
Kl — > 7r°yU + /i~ amplitude, with the correct chiral dynamics of the Kl — > 37r vertex and 
the missing kaon-loop contribution. In order to check the stability of this result, we also 
performed a detailed study of possible higher-order contributions, taking into account 
the recent experimental information on the Kl —>■ 7r°77 decay distribution. The main 
outcome of this analysis is a reliable theoretical estimate of the ratio 



Despite numerator and denominator receive, separately, large corrections beyond 0(p 4 ), 
i? 77 turns out to be a rather stable quantity, with a theoretical error conservatively esti- 
mated to be around 30%. 

Combining the analysis of the CPC amplitude with an updated estimate of the CPV 
terms, we finally arrive to the prediction B(Kl — > vr°/i + /i~)sM ~ 1-5 X 10~ n for the 
total branching ratio. This turns out to be composed by almost equal amounts of CPC, 
indirect-CPV and direct-CPV (including the interference) terms. The present error of 
this estimate is around 30%, but it is largely dominated by the parametric uncertainty 
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due to the B(Ks — > 7r°£ + £~) measurements. The irreducible theoretical uncertainty is 
around 10% and can be further decreased by suitable kinematical cuts. We thus conclude 
that the Kl — > 7r° / u + / u _ channel is a very promising mode for future precision studies 
of CP violation in the AS = 1 sector, with complementary virtues and disadvantages 
with respect to the other golden modes, namely Kl — > -n^vv and Kl — > TT°e + e~. As 
an illustrative example of the discovery potential of the two Kl — > tt°£ + £~ channels, we 
conclude our analysis showing that - even with the present parametric uncertainties - 
these two modes would yield a clear non-standard signature for the new-physics scenario 
recently proposed in Ref. [T2] . 

The paper is organized as follows: in Section 2 we discuss the kinematics and the 
general amplitude decomposition of Kl — > ir°C + Cr decays. In Section 3 we update the 
prediction of the CP- violating rate for the muon mode. Section 4 contains the main results 
of this work, namely the analytical and numerical study of the K L — > tt°{£ + £~)j=o CP- 
conserving amplitude. The phenomenological analysis of Kl — > 7r°/i + /i~ decays, within 
and beyond the SM, is presented in Section 5. The results are summarized in the Con- 
clusions. 



2 Amplitude decomposition and differential distribu- 
tions 

The most general decomposition of the Kl — > 7r°£ + £~ amplitude involves four independent 
form factors (S, P, V, and A) [13J: 
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m e u(p)v(p')S{z,y) + m t u(p)^ b v(p ! )P(z,y) 
+ (pk + p 7r ) fi u(p)^v(p')V(z,y) + (p K + P*)»u{p)^sv(p')A(z,y) 
The two independent kinematical variables can be chosen as 

(p + p') 2 - 2 p K -(p-p') 
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The variable y coincides with cos^, where 0£ is the angle between Kl and i spatial 
momenta in the dilepton center-of-mass frame, and the two independent ranges are 



4r| < z< (1-rv) 2 , 



T < y = cosOe < 1 



(5) 
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With these conventions, the general expression of the differential decay rate is 



rePt{z)z\S{z,y)\ + rfz\P{z, 



tfadj/ 2 12 7T 5 

+ fi(z)(l - %y 2 )\V(z, y)\ 2 + [#(*)(! - /^y 2 ) + 4r?(2 + 2r* - *)] \A(z, y)\ 2 

+ Ar 2 (l - rl)Re [A(z, y)*P(z, y)} + Ar 2 ^(z)fo(z)yRe [V{z, y)*S(z, y)\ } 



(6) 



The y variable is CP-odd and leads to a strong kinematical suppression. For this 
reason, it is convenient to expand the four form factors in a multipole series: 

F(z, y) = F (z) +yF 1 (z) + ... F = S,P,V,A. (7) 

The CP properties of the amplitudes are the following: the odd multipoles of S and the 
even multipoles of P, V and A, vanish in the limit of exact CP invariance. All short- 
distance contributions to the higher multipoles are negligible, since they cannot be gen- 
erated by dimension-six operators. To a good approximation, we can also neglect all the 
long-distance contributions generated by two-photon amplitudes in the CPV multipoles. 
As a result, we are left with five potentially interesting terms: 

A ,P : direct CPV amplitudes of short- distance origin; 

Vo : CPV amplitude which receives both short-distance (direct CPV component) and 
long-distance (indirect CPV component) contributions [Til I15j : 

So, Vi : CPC amplitude induced by the long-distance transition K L — > 7r°(77)j =0j2 — > 
tt°£ + £~; with the J PC = ++ two-photon state contributing only to So- 1 

As mentioned in the Introduction, these amplitudes turn out to be of comparable 
magnitude despite they are generated by very different dynamical mechanisms. Given 
the normalization in Eq. (J2J), the CPV short- distance multipoles are of 0(lm{V t *Vtd)/e 2 ~ 
1(T 3 ) within the SM, Vo ind is of 0(e ~ 1(T 3 ) and the two-photon CPC terms are of 
0(a/(47r) ~ 10~ 3 ). We stress that these five multipoles are the only interesting terms 
also beyond the SM, and that only the short- distance terms (A , P , V^ dir ) could receive 
sizable non-standard contributions. 

The differential z distribution taking into account only the five potentially interesting 
terms can be written as 

dT _ dT C pc _|_ ^r C py 

dz dz dz ' 



1 Note that this decomposition of the amplitude in terms of Dirac structures and y multipoles, which 
is quite useful to identify the dynamical origin of the various terms and their relative weight in the decay 
distribution, does not correspond to a projection on \£ + £~) states of definite angular momentum. Thus 
the So term receives both contributions from J PC = ++ and J PC = 2 ++ while, by construction, V\ 
receives only contributions from J PC = 2 ++ . 
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where 



dT CPC _ a*GWMz)fr(z) \ rimASQ{z)? + 2 ^ A + 6rf 



^ 2 n 7T 5 [ 15 ^V~; ^ z 

+^r|Ar(«)A(«)Re[5oWVi(«)l} , (9) 
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+4r 2 (l-rl)Re{A ( Z yP (z)] 

The asymmetric integration on y give rise to the following CPV distribution 
gL4fr /■ rf 2 r . a 2 G 2 F m 5 K (3 2 {z)(3 e {z) 



dz 



I dy -j-jz sgn(y) 
J dzdy 
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With a proper normalization, Ape can be identified with the forward-backward or the 
energy asymmetry of the two leptons. 



3 The CP-violating rate 

An updated discussion of the CP-violating amplitude for the electron case {I = e) can 
be found in Ref. |2|, to which we refer for more details. Most of the ingredients of this 
analysis can easily be transferred to the muon case. Following the notation of Ref. |2*lll5j. 
the dominant indirect- CPV amplitude leads to 

n±e(a s + b s z), (12) 

where the dz refers to the relative sign between this long-distance term and the short- 
distance (direct-CPV) component of the vector amplitude. 2 The latter is given by 

VS s '{z)=i^^U(z)lmX t , (13) 

where, as usual, X t denotes the CKM combination V t * s V t d, f±(z) are the K — > n form 
factors of the sj^d current, and yjv,A are the Wilson coefficients of the leading flavour- 
changing neutral-current (FCNC) operators ^Hj- Similarly, the axial-current operator 

2 The separation of V lnd and V dlr is phase-convention dependent. We adopt the standard CKM phase 
convention where, to a good approximation, we can set e = |e|e l7r / 4 . 



Vt d (z) = ±e 



a s + b s z + 
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leads to the clean short- distance terms 



A (z) = i^Af + (z)Im\ t , (14) 
P Q (z) = -i^f4z)lm\ t . (15) 

Note that, contrary to the electron case, here the pseudoscalar amplitude is not helicity 
suppressed (see Eq. (fTUj) ). 

Thanks to the recent NA48 results, 

B{K S -> ir°e + e~) mee>165 MeV = (3.0iJ;ij ± 0.2) x 10" 9 , (16) 
B{K S -> ttVV") = (2.9t\i ± 0-2) x 10~ 9 ^ , (17) 

the large theoretical uncertainty about the size of the indirect-CPV amplitude has been 
strongly reduced. Under the natural and mild assumption that the form factor of the 
K s — » n°£ + £~ amplitude is equal to f+(z), i.e. assuming (bs/ag)z = f + (z) — l = 0.39,2 [15J, 
the combination of these two measurements leads to 

\a s \ = 1.2 ± 0.2 . (18) 

Within the SM, the main parametric uncertainty on the direct-CPV amplitude is due to 
lm\ t . According to recent global CKM fits [T7j : 

ImA t = lm(V t * s V td ) (1.36 ± 0.12) x 10~ 4 . (19) 

Keeping the explicit dependence on these two main input values, and combining all the 
CPV contributions, we obtain the following phenomenological expression: 



B(K L -> ttV/Ocpv = 10" i2 x 



3.7|a 5 | 2 ±1.6|a s | — j ) +1.0 



ImA+ \ / ImA, x 



lO" 4 / ' VlO" 4 



• (20) 



The numerical coefficients have been obtained setting y 7 A = —(0.68 ± 0.03)a(M^) and 
yrv — (0.73±0.04)a(Mz), corresponding to m t {m t ) = 167±5 GeV and a renormalization 
scale between 0.8 and 1.2 GeV ^Hj. Using these inputs, the error on each of the three 
coefficients in Eq. (120 jl is below 10%. The result depends very mildly on /_, which 
has been set to f-(z) = —0.14; however, the relative weight of the three terms in (f2T)|) 
is different with respect to the electron case. This happens because of the terms 
in Eq. (jlOj). which enhance the axial-current direct-CPV contribution. As discussed in 
Ref . [2j , the preferred theoretical choice for the sign of the interference term is the positive 
sign. The same conclusion has recently been reached also by Friot, Greynat and de 
Rafael, who reinforced this theoretical argument with a detailed dynamical analysis of 
K — > 7T7 — > tt£ + £~ amplitudes in the large Nc limit 
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4 The CP-conserving rate 



As discussed in the Introduction, the leading CHPT contribution to the CP-conserving 
Kl — > 7T°£ + £~ amplitude is generated by two-loop diagrams of 0(e 4 p 4 ). The general 
structure of these diagrams is indicated in Fig. Due to the absence of local contri- 
butions, the complete two-loop amplitude of 0(e 4 p 4 ) is finite and does not depend on 
unknown coefficients. 

4.1 General decomposition of the loop amplitude 

From the general structure in Fig. [H we can decompose the CPC amplitude as 

A(K L - n°e+ncpc = ^ / 705 A i>' -4>P+ 90 ?') , (21) 

where 

L^(p,p>;q>)=u(p) 
and AOj" is defined by the Kl — > 7r°77 amplitude with off-shell photons: 

A[K L -> 7r°7(e 1 ,t)7(e 2 ,A;)] = 4£(t, . (23) 

A strong simplification arises by the observation that the leptonic tensor obeys the Ward- 
identity relations 

t v L» v = hpIS" = . (24) 

This implies that we can neglect in all the manifestly off-shell terms of the type tH v , 
kH u , k^k v . The non-trivial Lorentz structure of AftjZ is then identical to the one of the 
on-shell amplitude O ED] and its calculation proceeds in a very similar way. 

As shown in Ref. [Ej, employing a suitable basis for the pseudoscalar meson fields, it 
is possible to reduce the number of independent diagrams contributing to the Kl vr°77 
amplitude to the four terms in Fig. |21 We note here that this structure can be further 
simplified with a suitable decomposition of the K — > tt°P + P~ amplitude, which acts as 
an effective vertex in the diagrams EtHSt 1 - At 0(p 2 ) we can write: 

A[K -> n P + {p + )P-{p-)] = i [ao + a lZ + a 2 (p 2 + + p 2 _ - 2m%)] , (25) 

where z is defined as in Eq. © and the Oj do not depend on the meson momenta. At this 
order gauge invariance forces the effective vertices of the diagrams |2h-|2li to be completely 
determined by Eq. (|25jl via the minimal substitution. It is then easy to realize that the 
complete sum of the four terms in Fig. 121 is equivalent to the sum of EM2b performed 
without the off-shell term We are therefore left with the calculation of only two 
irreducible diagrams, with an effective K — > tt°P + P~ vertex which does not depend on 
the loop momenta: 

M% p {z) = a + a x z . (26) 
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Figure 1: General decomposition of the K L — > 7r°£ + £ amplitude. 




o a) 




o b) 
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Figure 2: Lowest-order diagrams contributing to the Kl 
of Ref. P2]. 



7r°77 amplitude in the basis 



This re-organization of the calculation leads to several advantages: i) a reduction of the 
number of relevant diagrams; ii) a manifest link with the Ks — > £ + £~ calculation of 
Ref. [TTj; hi) the possibility to include a well-defined class of higher-order contributions 
(see e.g. Ref. |2*T]). As we shall show later on, the latter step is achieved in the case of 
the dominant pion loops by the replacement of M^(z) with the K — ► 7r°7r + 7r _ amplitude 
determined by experiments. 

The explicit integration over the P-meson loop in AOZ, namely 



A%(p'-q',p + q 



-e z M y ; 

-7T U 



(2) 



pp\ 



d d q {2q - kY(2q + tf - g^ u {q 2 - mp) 
{2u) d [{q + t) 2 - m 2 P }[(q - k) 2 - m 2 P ][q 2 - mj 
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leads to 



A{K L -> tt 1 



0/ 



dp) 

>CPC 



x I dx I dy 
Jo Jo 



-—M {2) (z) 

d d q' u{p)^ v {— rri()^v{p') 
{2ir) d [{qiy - m 2 ]{q> + p) 2 {q> - p') 



x 



Axy -2xy(zm 2 K - 1 2 - k 2 ) + y(2y - l)t 2 + x(2x - l)k 2 



A(x,y) 



A(x,y) 



(27) 



where 



A(x, y) = m 2 P + ie — xy(zm 2 K ) + t 2 (— 1 + x + y) + k 2 {— 1 + x + y) . 



(28) 



4.2 Complete analytical expression and numerical integration 

The second loop integration is definitely more involved, but can still be expressed in a 
compact way in terms of four Feynman parameters. Using the equation of motion on the 
external lepton fields we can write 



A{K L -> 7T°£ H 



) {P) 
Jcpc 



cr 
2^ 



(2) , ,m t u{p)v{jf) ( rf r 



■M% P {z) 



zm 



K 



z 2; 



(29) 



/»1 pl—x rl rl 

1 (a, b) = dx ydy I dw I du 
Jo Jo Jo Jo 



where 



a + f3u — ie 



x \(2x + 2y-l)[l + u(l-w)]- X } , (30) 

1 + W) + IE 



a — w [b + y (y — 1)] , 7 = — l + w- 



wx 



(l-x-y)(x + y) 



P = wy[l ■ — (1 - x - y) + a (1 - w) (x + y) (1 - x - y) . (31) 

V x + yj 

The notation of Eqs. (|3*0j) and (|3*T|) has been chosen in order to show the perfect analogy 
and the agreement with the results of Ref. jTTj, in the appropriate kinematical limit. As 
discussed by Ecker and Pich, the integral (|3Uj) has three possible absorptive cuts: the 
77 cut, and the PP and PP7 cuts (both open only for b < 1/4). The real contribution 
associated to the double PP-77 cut turns out to be vanishing. 
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Figure 3: Behavior of l(m 2 /M 2 , m 2 /M 2 ) as a function of M (in MeV). 



These analytic properties become manifest performing the integration over the u vari- 
able in Eq. ([50]). This leads to 



Re Z(a, b) 



1 fl—X fl 

dx / ydy / dw 
a Jo Jo 

2x + 2y-l 
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log 
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Im Z(a, b) 
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Jo Jo Jo 



2x + 2y - 1 / a 

1 (1 — w 

2P V P 



\a + p\ 



a 4-/3 



4- lim 



a; (ck7 — /3) 1 / |a + (3\ 



H , |l+_7j _ x 
a I4-7 



(32) 



I a I 
a 



(33) 



The above representation, which holds independently of the possible signs of a, (3 and 
7, takes into account all the discontinuity prescriptions and is suitable for numerical 
integration. This last step is actually quite delicate because, though their combination 
is finite, both the PP and the PP7 cuts are infrared divergent (the electron case is 
especially difficult since collinear singularities are much stronger). This means that large 
cancellations occur between the various terms of X(a, b). We used a numerical integration 
routine based on VEGAS [22]; implemented in Mathematica [2B] on which we relied for 
preventing loss of numerical significance. The results are shown in Figs. OHI] for P = 71 
(where M < should be understood as i|M|). The behavior for P = K is entirely 
determined by the 77 cut and it is rather smooth (as can be guessed from Figs. OHH below 
2mA 



10 




Figure 4: Behavior of l(ml/M 2 , m 2 /M 2 ) as a function of M (in MeV). 

The values at z — 1, relevant for K$ — > 7r + 7r~ — ► 77 — > are 
T(r e V*) = 2.11-40.41z (±0.15) , 

J (rj,7^) = -2.821 + 1.2162 (±0.001), 1 ' 

which are compatible with the figures quoted in Ref. [TT] . To further check our results, we 
have taken advantage of the fact that A(J° (s) — > 7r°£ + £~) CPC satisfies an unsubtracted 
dispersion relation. In particular, the reduced amplitude X{M 2 ) = J(mj/M 2 ,mp/M 2 ) 
obeys the relation 

te EMjrjL,^. (35 ) 

M 2 it J s-M 2 s v 1 

After having computed the imaginary part over a large range of s from Eq. (jHSJ), we ob- 
tained a very good agreement between the real part values found by numerical integration 
of (|35|) or directly from ([32)1. for both M 2 > and M 2 < 0. This check is particularly 
significant in the region M 2 < 0, where the integration of (|3*2j) is stable, being singularity- 
free, while (135)) remains very sensitive to the cancellations among PP and the PPj cuts 
in ImX(s). 

4.3 Total and differential rate, and inclusion of higher-order 
terms 

Replacing the 0(p 2 ) amplitude Mq pp (z) in Eq. with the following parameterization 

M% p {z) - G 8 m 2 K a p (z) , (36) 
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Constant 


Chiral 


Dalitz 


a n {z) 
a K {z) 


-0.362 



2 

z-rl-l 


-0.46 + 2.44,2 - 0.95z 2 



Table 1: Phenomenological functions a p (z) used in the — > tx°P + P vertices. 







No KK 


Chiral KK 






Constant txtx 


Chiral txtx 


Dalitz txtx 


Dalitz txtx 




£>CPC 


3.48 x 10~ 12 


2.26 x 10- 12 


2.62 x 10- 12 


2.81 x 10~ 12 






4.15 x 10" 6 


3.80 x 10~ 6 


3.58 x 10" 6 


3.66 x HT 6 




Rabs 


81% 


79% 


77% 


84% 


e + e~ 


£>CPC 


2.28 x 10" 14 


1.48 x 10" 14 


1.75 x 10" 14 


1.80 x 10" 14 






2.73 x 10~ 8 


2.47 x HT 8 


2.39 x 10~ 8 


2.34 x HT 8 




Rabs 


27% 


28% 


30% 


28% 



Table 2: Branching ratio, ratio _R 77 and absorptive component, R a b s = r^pc/rcpc, for 
the two Kl — > tx°£ + £~ modes, as obtained with different phenomenological expressions of 
a*> K (z). 



where G$ 

dTcpc 
dz 



9.1 x 10~ 6 GeV 2 , the CPC differential rate becomes 
G\ra\a A r 2 



5127T 7 



z z 



r 2 r? 



,K 



z z 



(37) 



We have evaluated this expression using several parameterizations of a p (z), reported in 
Tabled The constant one corresponds to the first approximation of Ref. ^U]. The chi- 
ral forms, deduced from the O (p 2 ) amplitudes for K — > tx°tx + tx~ and K — > tx°K + K~, 
are the expressions which yield to the complete 0(p 4 ) calculation of the CPC rate. Fi- 
nally, the Dalitz expression is obtained from the experimental fit of the Kl — > tx°tx + tx~ 
amplitude |24ll25]: 



A{K L ^ 

a(z) 



TX TX 



V) = Gs m 2 K a(z) + O {[p K ■ (p + - p_)] 2 } 



0.38 + 0.13F - 0.0059F C 



) 



Yn= — 



(3? 



neglecting the small O {[pk • (p+ ~ P~)} 2 } terms. The resulting branching ratios for the 
muon and electron modes are shown in Table 121 The differential rate dT /dz with or 
without the kaon loops and with the Dalitz parametrization for a n is shown in Fig. As 
can be seen, the impact of the charged kaon loops is quite small. 

Of special interest is the ratio _R 77 , defined in Eq. (JTJ), with the two widths computed 
in terms of the same a p (z) vertex. Concerning the two-photon mode, this means j2I 



T (K L -> 7T 7 7 ) 



Glm 5 K a 2 



(4tt)' 



dz f3 n (z) |a 7r (z)X 77 (z/rl) — a (z)T^(z) | 



(39) 
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dr CP c 5 

dz 4 

3 
2 
1 

" 0.2 0.3 0.4 0.5 * 

Figure 5: Differential rate for the CPC contribution, in arbitrary units, including (plain) 
and excluding (dashed) the K + K~ loop. 

with the standard loop function Z 77 given in Ref . [TH] . 

The phenomenological functions a p (z) act as a modulation of both the £ + £~ and 77 
spectra, and the four cases in Tabled span a large range of behaviors (this is why we kept 
the constant a p , even if this choice is clearly in conflict with chiral symmetry and the 
experimental Kl — > tt 77 spectrum). Interestingly, the dynamics of the processes is such 
that Kl —>■ 7r°77 and Kl — > n°£ + £~ amplitudes react in a very similar way to changes 
in the distribution of the invariant mass of the electromagnetic system (i.e. the 77 or 
£ + £r pairs). As a result, despite the individual variations of the two widths, the ratio i? 77 
remains very stable. Note that this conclusion is not induced by a complete dominance of 
the 27-cut in the Kl — > 7r°£ + £~ amplitude: tot and 7T7T7 cuts have sizable contributions, 
and the stability of _R 77 holds also for the Kl — > 7r°e + e _ amplitude, which is dominated 
by the dispersive part. 

As is well-known, none of the parameterizations in Table 1 for a p is able to reproduce 
both rate and spectrum of the K L — > 7r°77 decay, and some additional O (p 6 ) genuine 
local terms need to be included (see e.g. Refs. EH] for a recent discussion). As shown 
in Ref. [2], recent data point toward O (p 6 ) counterterms whose bulk values can be inter- 
preted in terms of vector-meson resonances. The overall effect of these extra contributions 
is an enhancement of the Kl —>■ 7r°77 rate, with a spectrum which remains very similar 
to what we have called the Dalitz structure. Unfortunately, this information cannot be 
directly translated into a complete O (p 6 ) analysis of the Kl — > tt°£ + £~ amplitude. How- 
ever, given the stability of _R 77 in Table El we expect that _R 77 will not change substantially 
also with the inclusion of these extra O (p 6 ) terms. This expectation is further reinforced 
by the observation that the Kl — > 7r°77 spectrum does not show significant evidences 
of two-photon states with J 7^ [2], and by the large size of i? a b s in Table El Since we 
cannot fully determine the structure of the Kl — > n £ + £~ amplitude at 0(p 6 ), we shall 
attribute to our final estimate of i? 77 a conservative 30% error, dictated by naive chiral 
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power counting. For the interesting muon mode, this means 



J R 77 (£ = /i ) = (3.7 ±1.7) x 10' 



(40) 



where the central value has been obtained from the last column in Tabled (Dalitz expres- 
sion for the pion loops and inclusion of the kaon loops). 

In order to obtain the predictions of the CPC Kl — > 7t°£ + £~ widths, the estimates 
of -R 77 need to be combined with the experimental results on B (K L — > 7r°77). The two 
most recent determinations 



B (K L -> 7r 77 ) 



cxp 



;l.36 ± 0.03 stat 
;i.68±0.07 stat ±0.08 



± 0.03 syst ± 0.03 norm ) x 10' 



syst, 



x 10 



-6 



5J 
6J 



(41) 



leads to the average B (Kl — 
we finally obtain 

B(K L 
B(K L 



7T°77) 



cxp 



^1.42 ±0. 13) x 10 6 . Using this averaged result, 



o + -\0++ 

7r e e J CPC 

+ -\0++ 
7T» ) CPC 



(3.4 ± 1.0) x 10~ 14 
(5.2 ± 1.6) x 10~ 12 



(42) 



Given the strong helicity suppression, the result in Eq. (J42j) is not the dominant CPC 
contribution for the electron mode. As shown in Ref. [2J, the constraints on the Kl — ► 
7r°77 diphoton spectrum (at low invariant mass) leave room for a 2 ++ CPC branching 
ratio, in the electron channel, around or below the 10 -12 level. Because of the phase-space 
and angular-momentum suppression, this translates into a 10~ 14 contribution from 2 ++ 
in the muon mode, which can thus be safely neglected. 



5 Phenomenological analysis 

The total branching ratios for the two modes can be written as 

B(K L -> *W ) = (c^ ± CL (5^) + Ci x (^j 2 + CS PC j x 10- 12 , (43) 

with the following sets of coefficients (obtained by combining the results of the previous 
two sections and the analysis of Ref. j2] for Kl — > rr°e + e~): 



C* ix = (15.7±0.3)|a 5 | 2 , C^ ix = (3.7 ± 0.1)|a s | 2 
Cf nt = (6.2±0.3)M, Cf nt = (1.6±0.1)H, 



Q ir = 2.4 ±0.2, = 1.0 ±0.1, 

^cpc ~ j ^cpc = 5.2 =t 1-6 , 

-4\ 



(44) 



where \as\ = 1.2±0.2 and (ImAi/10 ) SM = 1.36±0.12. The Standard Model predictions 
are then 



B(K L - 7r e+e-) = (3.7±S;J) x lO" 11 



B(K L -> TT°fi + /2-) = (1.5 ± 0.3) x 10~ n 



(1.71^) x 10-") ] , 

(1.0 ±0.2) x 10' 11 ) ] , { > 
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dr 
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Z 
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0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Figure 6: Differential rate for Kl — > 7r°/i + / u _ , in arbitrary units, for destructive (left) 
and constructive (right) interference. In each plot the three curves correspond to CPC 
(dashed), CPV (dotted) and total contribution (full). 

for positive [negative] interference. As already mentioned, several theoretical arguments 
point towards the constructive interference scenario (21 CH] ■ We report the destructive 
solution only for completeness, since we cannot prove yet in a model-independent way 
that this solution is ruled out. 

The total differential rate and its CPC and CPV components are shown in Fig. 13 for 
the muon mode (the electron case is less interesting, being completely dominated by the 
CPV part). Although the differential CPC component is not predicted as accurately as 
the total one (by means of -R 77 ), it clearly emerges that the low fi + n~ invariant mass 
region is dominated by the CPV component. This conclusion can be considered as a 
model-independent statement, at least at a qualitative level: it follows from the naive 
distribution of |/i + /i~> final states with J CP = 1 ±± (CPV) or 0++ (CPC), and it is 
reinforced by the experimental evidence of the — > 7r°77 diphoton spectrum |5J Ej. 
For this reason, in future high-statistics experiments a cut of the high-z region could be 
used as a powerful tool to reduce and control the CPC contamination of the signal. At a 
quantitative level, employing the Dalitz shape of the spectrum, we find that cutting the 
events with m MAt > 2m n let us suppress about 90% of the CPC contribution, while the 
interesting CPV term is reduced only by 40% (i.e. the CPC contamination of the total 
rate drops to less than 10%). 

The different shape of CPC and CPV spectra implies a suppression of their interference 
in the CPV distributions, such as the forward-backward (FB) asymmetry. In particular, 
we find that the normalized FB asymmetry integrated over the full spectrum do not 
exceed the O(10%) level. Thus we do not consider this observable particularly promising 
for near-future (low-statistics) experiments. On the other hand, it is clear that in a 
long term (high-statistics) perspective this observable, as well as the other asymmetries 
discussed in Ref. [27], would provide a useful tool to reduce the theoretical error on the 
long-distance components of the K L — > 7r°/i + /i~ amplitude. 
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Figure 7: Central values of the two B(Ki — > ir°£ + £ ) for various values of Im X t . Negative 
values stand for destructive interference (see text). 

The present experimental limits on the branching ratios of the two modes can be 
translated into bounds on lmX t , although we are clearly far from the level necessary to 
perform precision tests of the SM. We find 

B{K L -> rr°e + e-) < 2.8 x 1(T 10 j2H] -> |ImA t | < 1.3 x 1(T 3 , . , 

B (K L -> ttVaO < 3.8 x 10" 10 [2H] -> |Im A t | < 2.1 x lO" 3 , ^ 

where the results have been obtained without any assumption on the sign of the interfer- 
ence. It is worth to stress that the limit derived from the muon mode is quite close to the 
electron one. This happens because the phase-space suppression of the former is partially 
compensated by the enhancement of the helicity-suppressed terms in the (direct-CPV) 
axial-current contribution. 

5.1 New-physics sensitivity 

As can be deduced by the comparison of Eq. and Eq. ([45)1 . we still have to wait a few 
years before any experiment will be able to probe the two Kl vr°£ + £ _ modes at the SM 
rate. However, we recall that these two transitions are particularly sensitive to physics 
beyond the SM and their branching ratios could well be modified (and possibly enhanced) 
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Figure 8: 25%, 50% and 75% confidence- level regions for the Standard Model (assuming 
positive interference) and the enhanced-electroweak-penguin model of Ref. taking 
into account all the present uncertainties. 

with respect to their SM expectations (see e.g. Refs. [3U] and references therein). In such 
a case, the combined information on the two rates would provide a very powerful tool to 
identify a deviation from the SM, and also to distinguish various new-physics scenarios. 

To illustrate the combined discovery potential of the two modes, it is useful to draw 
the curve in the B(K L — > 7r°/i + AO~ B(K L — > 7r°e + e~) plane obtained by a variation of 
ImAt, as shown in Fig. The distance between the positive and negative ImAt branches 
is generated by the different weights of the CPV contributions for the two modes, which in 
turn arise from the helicity-suppressed terms proportional to Pq(z) and A (z) in Eq. (JTUJl . 
In the limit where we neglect all errors but the one on ImAt, as done in Fig. [7| for 
illustration, the SM corresponds to one small segment (or two considering also the negative 
interference) of this curve. Any other point along the curve corresponds to non-standard 
scenarios where the new-physics effect can be re-absorbed into a re-definition of ImAt. 
Finally, the region outside the curve corresponds to non-standard scenarios with different 
weights of vector and axial-vector contributions. 

To discuss a concrete example, we shall consider in particular the scenario with en- 
hanced electroweak penguins recently discussed in Ref. [T2]. In this framework new 
physics does not modify the value of ImAt, but leads to a modification of the short- 
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distance Wilson coefficients yrv,A = VivaI a ^Mz)- In particular, their values are changed 
to (yiAiViv) — (—3.2,0.9), whereas the SM case corresponds to (—0.68,0.73). Given the 
dependence of Cf nt dir from yrv,A 

C[ nt = 8.91 y 7V \a s \, = 2.67 (f 7V + y 2 7A ) , 

tf& = 2.12 y 7V \a s \, = 0.63 (y 2 7V + y 7 2 A ) + 0.85 y 2 7A , ^ 

we can easily translate this information into the corresponding predictions for the two 
branching ratios: 

B (K L - 7r°e + e-)^ P = (9.0 ± 1.6) x lO" 11 , 
B(K L - tt W)ewp = (4-3 ± 0.7) x lO" 11 . 

Both decay widths turn out to be dominated by the direct CPV component, with a 
small (20%-30%) dependence from the interference sign. Since the helicity suppressed 
effects are proportional to the enhanced y 7 A, the relative deviation with respect to the 
SM prediction is larger in the muon case. As shown in Fig. |HJ this scenario could clearly 
be distinguished from the SM case even taking into account all the present parametric 
uncertainties. Note that, as a consequence of the non-standard axial/vector ratio, the 
central point of this new-physics scenario is outside the Im X t curve. 



6 Conclusions 

We have presented a new comprehensive study of the rare FCNC process K L — > 7r°// + /-*~- 
Our main task has been a systematic analysis of the two-photon CP-conserving contri- 
bution, which represents the most delicate ingredient necessary to obtain a reliable SM 
prediction for the total decay rate. By means of an explicit two-loop calculation within 
CHPT, we have shown that this contribution can be predicted with reasonable accuracy 
in terms of the Kl — > 7r°77 width. Taking into account the experimental data on the 
latter jSJ |Hj and the recent observation of the Kg — > ti £ + £~ transitions [3 E] , we have 
been able to present the first complete SM prediction for the total branching ratio: 

B(K L -> ttV + /0 SM = (1-5 ± 0.3) x 10- 11 . (48) 

This is composed by approximately equal CP-conserving, indirect-CP- violating and direct- 
CP- violating components (including the interference). The sizable uncertainty which af- 
fects this prediction is a parametric error which reflects the present poor experimental 
knowledge of the Kg — > tt°£ + £~ rates. The irreducible theoretical error, due to the long- 
distance CP-conserving amplitude, is only at the 10% level and, as we have shown, can 
be further reduced by suitable kinematical cuts. 

Beside the good theoretical control of the decay rate, this channel is particularly 
interesting because of a different short-distance sensitivity compared to the other golden 
rare Kl decays, namely Kl — > 7r°e + e _ and Kl — ■> vr°z/z/. As we have shown, the helicity- 
suppressed terms distinguish the relative weights of vector- and axial-current contributions 
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in Kl — > TT°e + e~ and Kl — * TT°fi + fi~. As a result, the combination of the two decay widths 
could provide a very powerful tool to falsify the SM and also to distinguish among different 
new-physics models. 

In summary, the K L — > 7r°fi + fi~ decay represents a very interesting candidate for future 
precision studies of CP violation and new physics in AS" = 1 transitions. This mode has 
complementary virtues and disadvantages compared to the other rare K decays, and it 
should be seriously considered in view of future high-statistics experiments in the kaon 
sector. 
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